Noise in an ac biased junction. 
Non-stationary Aharonov-Bohm effect. 



G.B. LesovikB 

Universitat zu Koln, 
Institut fur Theoretische Physik, 
Ziilpicher Str. 77, D-5000 Koln 41, Germany 

L.S.LevitovS 

Massachusetts Institute of Technology, 
12-112, Department of Physics, 
77 Massachusetts Ave., Cambridge, MA 02139 



Abstract 

We study excess noise in a quantum conductor in the presence of constant 
voltage and alternating external field. Due to a two particle interference effect 
caused by Fermi correlations the noise is sensitive to the phase of the time 
dependent transmission amplitude. We compute spectral density and show 
that at T = the noise has singular dependence on the dc voltage V and the 
ac frequency f2 with cusplike singularities at integer eV/fiQ,. For a metallic 
loop with an alternating flux the phase sensitivity leads to an oscillating 
dependence of the strengths of the cusps on the flux amplitude. 

PACS numbers: 72.10.Bg, 73.50.Fq, 73.50.Td 
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There is a variety of phenomena related with the quantum coherence of transport in 
small conductors 0: weak localization, Aharonov-Bohm effect with the flux quantum hc/2e, 
universal conductance fluctuations, etc. Each of these effects can also be seen in the spectrum 
of noise, equilibrium or non-equilibrium. The equilibrium noise is simply proportional to 
conductance according to the fluctuation-dissipation theorem. The non-equilibrium noise in 
coherent conductors is expressed through eigenvalues of the scattering matrix fHHH, and, 
therefore, is also related with the conductance, though in a less trivial way. For that reason 
all the coherence phenomena are present in the non-equilibrium noise as well. However, for a 
better understanding of transport in small conductors it is interesting to analyze the converse 
line of thinking and to look for coherence effects that are present in the noise but are absent in 
the conductance. Such effects, if they exist, are genuinely many-particle ||, otherwise they 
would show up in the conductance. As long as we are talking about non-interacting fermions 
it is only statistics that can produce such coherence. A purpose of this letter is to describe 
an effect caused by two-particle statistical correlations that leads to phase sensitivity of a 
two-particle observable, i.e., of electric noise, but does not affect one-particle observables, 
e.g., conductance. The phase sensitivity manifests itself in an oscillating dependence on the 
amplitude of an ac flux, in many aspects similar to the A-B effect. However, it will occur 
in a single-connected conducting loop, i.e., in the geometry where the normal A-B effect is 
absent. 

Let us specify which coherence effects we are going to study. In simple words, when an 
electron is scattered inside a conductor its wavepacket splits into two portions, forward and 
backward, presenting a choice to the electron to be either transmitted or reflected with the 
probabilities D and 1 — D. Part of this picture of the wavepacket splitting, involving the 
relation of D with the conductance ||[7j and of D(l — D) with the noise 0J^,|,|J is well un- 
derstood. However, there is another part, quite unusual, related with the behavior of current 
fluctuations in the time domain. Recently, we studied the distribution of the charge trans- 
mitted through a resistor during fixed interval of time M. We found that the distribution is 
very close to the binomial, which means that the attempts to have electron transmitted are 
highly correlated in time. (Were the sequence of the attempts perfectly periodic the distri- 
bution would be exactly binomial.) The origin of the correlation is the Pauli principle that 
forbids passing of electrons through the resistor simultaneously. The attempts follow almost 
periodically, spaced by the interval h/eV. Because the periodicity is not perfect it does not 
affect the average current, but shows up in its second moment, i.e., noise, leading |L0| at 
zero temperature to a sharp edge of the spectral density of excess noise near u = eV/h: 

= 2^D(1 — D)Ti{ujq — \uj\) for |u;| < u®, otherwise. (Excess noise is the difference of the 
actual noise and the equilibrium S u = 2^-Dhu coth(huj/2T).) The corresponding current- 
current correlation in time is {(J{t)j(t + r))) = ^-D(l — D)sm 2 [ujqt / 2) / r 2 , oscillating with 
the period 2ti/lj and decaying. 

Having realized that the frequency eV/h is characteristic for the time correlation of 
the attempts one has to think of a simple experimental situation were the presence of this 
frequency could be studied. It is natural to consider a system driven both dc and ac, and 
to look for the effects of commensurability of Q and eV/h, where Q is the frequency of the 
ac bias and V is the dc voltage. In this letter we study such a system and demonstrate that 
due to the ac bias the singularity at u = u Q can be shifted down to zero frequency thus 
making it easier to observe. Below, we compute the noise in a model resistor in the presence 
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of combined dc-ac bias and find that the low frequency noise power Sq has singularities 
at eV = nhfl, when the "internal" frequency of the problem eV/h is a multiple of the 
external frequency Q. We find that 8S /dV is a stepwise function of V that rises in positive 
steps at V n = nktt/e. Another interesting observation is that the heights of the steps of 
dSo/dV are phase sensitive, i.e., they depend on the phase of the transmission amplitude 
in an oscillating way resembling A-B effect. The phase sensitivity of the noise should be 
opposed to the pure dc situation where only the probabilities of transmission and reflection 
enter the expression for the noise, which makes the noise power insensitive to the phase 
picked by the wavefunction across the system. In the simplest situation when the ac bias 
is supplied by alternating flux threading the current loop, $(£) = $ a sin(fit), the heights of 
the steps in dSo/dV are proportional to the squares of the Bessel functions J^(27r<E> a /<I?o), 
where $o — hc/e. Let us note that we are not talking about the trivial effect of the e.m.f. 
—d^/cdt induced in the circuit by the alternating flux. The effect in the noise will persist 
in the quasistatic limit \d&/cdt\ <C V when the ac component of the current vanishes. 

Let us start with recalling general facts about scattering off an oscillating potential. We 
consider a model one dimensional system where electrons are scattered by alternating scalar 
and vector potentials U(x,t), A(x,t) localized in the interval [— d, d], U(x,t) = A(x,t) = 
for |x| > d. As a function of time they are periodic: U(x,t) = J2m=-oo U m (x) exp(— imflt), 
where U (x) is the static part of the potential, and the other harmonics U m (x), m ^ 
describe ac bias. (Expression for A(x, t) is similar.) The dc bias is expressed in the framework 
of the Landauer model as the difference of the population of the right and the left scattering 
states. An important difference is that in our case the states describe inelastic scattering 
because an electron can gain several quanta hQ while passing through the region [— d, d]. 
It will be useful to have the states expressed through the amplitudes of transmission and 
reflection: 

( e~ iEt+ikx + £ B Ltn e~ iEnt - iknX x < -d 

^L,fc(M) = j ^ ALne -iE n t+ik n x x>d , 



£ A R>n e~ iEnt - iknX x < -d 

*pR,k( X i = { 1-iEt-ikx + £ BRne -iE n t+ik n x x > d , (!) 
I n 

where the amplitudes ALm\ n ,Bi,m\ n are time independent. Here E = h 2 k 2 /2m, E n = 
E+nhQ, and k n are defined by f?k\ = 2mE n . The states (1) are solutions of the Schrodinger 
equation 

E^(x,t) = [^(-t^--^A(x,t)) 2 + U(x,t)}tP(x,t) . 

They can be used as a basis to study transport through the system the same way it is done 
for the static barrier. The amplitudes A L ( R ^ n , B L ( R ^ n satisfy unitarity relation, 

£ 5(E n - E' n ,)(A a , n (E)A a , tn ,(E') + B a , n (E)B a ,, n ,(E')) = 5{E - E')8 aa , 
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that one obtains by the standard reasoning about conservation of current. 

The operator of electric current, j(x,t) = —ieip + (x,t)Vip(x,t) is written in terms 
of second-quantized electrons, i>{x,t) = ipL(x,t) + ip R (x,t), ipL{x,i) = J2ki J L,k(x,t)ak, 
i/j R (x,t) = Y^k^Rtkixjtpk, where and bk are canonical Fermi operators corresponding 
to the states (1) coming out of the reservoirs, the left and the right respectively. It is 
straightforward to compute the mean value I(t) = (j(x,t)), where the brackets (...) stand 
for averaging with the density matrix p of the reservoirs. As usual, we assume absence of 
correlations in the reservoirs, p = Pl ® Pr, which physically means that after having been 
scattered into a reservoir electrons have enough time to relax to the equilibrium before they 
return. Below we assume equilibrium Fermi distributions pl )R = n(E — Ep ± eV/2). One 
obtains a generalized Landauer formula: 

oo 

I{t) = ]T J m exp(-Mt) , (2) 

m=— oo 

where it is straightforward to write the coefficients I m in terms of the scattering amplitudes 
Al(r),u and B L ( R ^ n . Expr.(2) describes both steady current and generation of harmonics in 
the presence of the ac bias. 

Quite similarly one can obtain an expression through A L ^ R ^ n , B L ( R ) n for the noise. 
Noise is related with the correlation function S(ti,t 2 ) = (([j(x, t\),j(x, £2)]+))- In the usual 
dc situation S(ti,t 2 ) = S(ti — t 2 ) and its Fourier transform gives spectral density = 
{(jwj-u)) of the noise. In our ac case the situation is somewhat more complex because the 
spectral density S u does not provide a complete description of the noise. Indeed, S(ti,t 2 ) 
will depend now separately on t\ and t 2 , not only on the difference t\ — t 2 . However, it 
satisfies S(ti,t 2 ) = S(ti + 2n/Q,t 2 + 2tt/Q) resulting from the periodicity of the ac bias. 
For the Fourier components j w it means that the average ({juju')) does not vanish whenever 
uj + uj' = mfl, where m is any integer. Thus, in addition one gets generalized spectral 
densities S^ )m = ((juijmSi-Jfli an integer parameter family of functions. Among them there 
is the 'ordinary' m = spectral density S u = S^fi = ((ju>j-u>)), the one easiest to access 
experimentally. In what follows we concentrate on it and do not study other S^m, m ^ 0. 

To compute the noise one has to average the product of two current operators over the 
distribution in the reservoirs. Evaluation of the average is similar to Refs. U0,||] , so we do 
not need to repeat it here. General expression simplifies quite substantially in the practically 
interesting limit of £/, the time of flight through the barrier U (x, t) being much shorter than 
2tt/Q and Ti/eV . The point is that h/tf defines the characteristic scale of energy dependence 
of the scattering amplitudes, so the condition tfQ 1, tfeV h enables one to neglect the 
energy dependence of A L ( R \ n , B L ^ >n in the interesting energy domain E F ±max[eV, HQ}. We 
also assume Ep 3> maxfeV, HQ], which allows to neglect the difference of k n /m and k/m, the 
velocities of scattered and incident states, and set k n /m = vp. It should be remarked that 
the physical picture we discuss below is not really dependent on any of these assumptions, 
they only make our expressions more compact. The more general case of arbitrary relation 
between h/tf, Ep, eV, hfl presents no difficulty. 

With the above assumptions made it becomes convenient to use Fourier transform of the 
amplitudes A and B. Let us define A a (t) = J2 n A a ,n exp(— inflt), a = L, R. Similarly we 
introduce B a (t), and rewrite Expr.(l) as 
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^L,k( X ^) 



e ikx + Bl( j + x / VF y-ikx x< _ d 
A L (t - x/v F )e ikx x > d 



. , _ j A R (t + x/v F )e ikx x < -d ( , 

V>R,kW, V ~ j e -ikx + Br( j. _ X / Vp ) e ikx x>d ■ W 

(To obtain (3) from (1) we substitute k n = k+nVl/v F in the phase shifts e lknX and then do the 
sum over n.) The amplitudes A L ^(t), B L ^(t) have clear meaning of the transmission and 
reflection amplitudes at given instant of time for a slowly varying potential. The retarded 
time t — \x\/v F in Expr.(3) accounts for the finite speed of motion after scattering. The 
unitarity relation now takes the form 

\A L(R) (t)\ 2 + \B L{R) (t)\ 2 = 1, A L (t)B R (t) + B L (t)A R (t) = . 

To clarify the character of the simplification thus achieved let us remark that with Expr. (3) 
the formula (2) for the current I(t) becomes just I(t) = 2^\A(t)\ 2 eV which means that the 
current 'adiabatically' follows time variation of the transparency of the barrier according 
to the Landauer formula. Now we shall compute noise and find that, unlike I(t), it is not 
reduced to anything trivially related with the static limit. Let us write the average of two 
currents ((jOOjfe))) = 

2£ £ e ^-W*-V[\A(t 1 )A(t 2 )\\n L (E')(l - n L {E)) + n R (E')(l - n R (E))) 

11 E,E> 



+B(t l )A(t l )A(t 2 )B(t 2 )n R (E')(l - n L (E)) + A^B^B^A^n^E!)^ " n R (E))]. 

To compute we have to do Fourier transform and substitute Fermi distributions 
til{r){E). Explicit calculation yields the result 

S u = — £ 2N o(^ - ^)l(l^| 2 )n| 2 + N^u.nn + eV)\(AB) n \ 2 , (4) 
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where 



Nq(x) = J (n(E -x) + n{E + x))(l - n{E))dE = xcoth(a;/2T) , 

a;sinh(a;/2T) — ysinh(y/2T) 



N 1 (x,y)=N (x + y) + N (x-y) = 



cosh(a;/2T) - cosh(?//2T) 



and (...)„ denotes Fourier components, e.g., (AB) n — ^ J A(t)B(t)e mQt dt. Expr. (4) de- 
scribes the noise as function of eV, f2, u and T. The behavior is simplest at T = when 
N (x) = \x\, Ni(x,y) = \x + y\ + \x — y\. Given by Expr. (4) as a weighted sum of terms 
like \nfl + eV ± ou\, \u — nfl\ the noise S u will then depend on V, Q, u in a piecewise linear 
way, changing from one slope to another when nfl + eV ± uj or uj — nfl equals 0. This 
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condition defines the locations where has singularities. They are cusps, sharp at T = 
and rounded on the scale T at T > 0. 

With the general Expr.(4) one can explore the noise in all possible limiting situations 
that one obtains for different combinations of eV, £1, u and T. Particularly interesting for 
us will be the case T = 0, uj = corresponding to the noise So = ((juj- u ))u-*o measured at 
low frequency. Let us remark here that setting u — means only that u is small compared 
to the parameters eV and Q that define the width of the frequency band of the excess noise. 
Such uj may still be much higher than the band width for other sources of noise, e.g., the 
1/ /. Let us concentrate on the dependence of So on V. It is a piecewise linear function 
which is easiest to characterize by its derivative, 

dSo/dV = — V \ n 6(eV - nhil) , (5) 

where A n = {(AlBr)^ 2 and 9(x) = 1 for x > 0, —1 otherwise. The function dSo/dV rises 
in positive steps at all V n = TiVLn/e (see Fig. |l|), the property that can be alternatively 
formulated as convexity of Sq(V) in V. 

The meaning of the singularities in Sq(V) was clarified recently in a study of the statistics 
of transmitted charge The generating function of the charge distribution was expressed 



through the single-particle scattering matrix, and it was found that the distribution arises 
from Bernoulli statistics (i.e., it is a generalized binomial distribution). The frequencies 
of attempts were given as function of V and Q. The probabilities of outcomes of a single 
attempt were found in terms of many-particle scattering amplitudes, and it was shown 
that they change at the thresholds V n = nhfl/e in a discontinuous way due to statistical 
correlation in the outgoing channels of the scattering. The discontinuity manifests itself in 
the second moment of the distribution that corresponds to the noise So(V) discussed above. 

There is an interesting and simple example where one can explicitly evaluate the heights 
of the steps. Let us consider a junction with ideal leads bent into a loop of length L 
(see inset of Fig. [I]) and placed into an external magnetic field varying with time. In this 
problem the junction is the only source of scattering. For simplicity let us assume that 
only one scattering channel is involved and that the junction is symmetric, Ai = A R = A, 
Bl = B R = B. The ac bias is supplied by the alternating flux of the magnetic field through 
the loop, $(t) = $ a sin(fit). Also let us suppose that the magnetic field is quasistatic, i.e., 
the time of flight through the system, tf = L/vp is much shorter than 27r/Q, that makes it 
possible to introduce the time dependent amplitudes A L ( R )(t), B L ( R ^(t) as it was discussed 
above. As is common, in such a situation the vector potential can be treated semiclassically, 
and one can write the wavefunction as if>(x,t) = exp(|| Jf^ A(x')dx')ipo(x, t), where x is 
the coordinate along the lead and ipoi^x^t) is found by solving the Schrodinger equation in 
the absence of the magnetic field. Thus all the dependence on the magnetic flux can be 
accumulated in the phase of the transmission amplitude, 

A m (t) = exp(±z$(t)/$ )A B R[L) (t) = B, 

where $o = hc/e is single electron flux quantum. Since |A(t)| 2 = D = const the cur- 
rent is time independent: / = ^-DV. According to Expr.(4) S u is written through the 
Fourier components of A L {t)B R {t) in this case given by the Bessel functions: (A L B R ) n = 
J n (27r$ a /$o)AB. Thus we find 
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2e 2 



[2N (u)D 2 + £ Ni(u, nQ + eV)D{\ - D)J 2 n {2^J%)}. 



(6) 



n 



The heights X n of the steps in dS /dV are then given by 



\ n = D(l-D)J 2 (2n$ a /%). 



(7) 



They oscillate as function of $ a / $o an d vanish at the nodes of Bessel functions. 

Exprs. (6) , (7) illustrate one important feature of the noise in the ac biased system, the 
sensitivity to the phase of the transmission amplitude A. By varying the amplitude $ a 
of the alternating flux one can make X n vanish separately for each harmonic nQ of the ac 
frequency. This should be compared with the case of the dc bias where the noise is expressed 
only through \A\ 2 and thus cannot be phase dependent. We call the oscillating dependence 
(6), (7) non- stationary Aharonov-Bohm effect. To compare it with the usual A-B effect let 
us recall that the latter is observed as an oscillation of the dc conductance under variation 
of flux in the situation when one has interference of transmission amplitudes corresponding 
to different classical trajectories of a quantum system, e.g., in a conductor with multiply 
connected leads forming one or several closed loops. The dc A-B effect cannot be observed 
in the single path geometry like Fig.l. Alternatively, the non-stationary A-B effect appears 
as a result of interference of the right and left scattering states travelling in the opposite 
directions along same path and having energies shifted by nQ. It is clear from our discussion 
that such interference does not contribute to the ac conductance but is important for the 
noise and, therefore, one obtains the non-stationary A-B effect in the noise even in the 
topologically trivial situation of Fig.l. 

One can derive a sum rule: 



where D(t) = \A(t)\ 2 . For A n given by Expr.(7) it follows from the definition of the Bessel 
functions. In the general case of Expr.(5) the sum rule is obtained by applying Plancherel's 
formula to Fourier components of A(t)B(t). The sum rule clarifies the relation of our problem 
with the previous calculation [@,[3],[|| of the noise in the pure dc case for which the result 
does not depend on the phase of A(t). When the limit is taken Q — > 0, V = const, the steps 
in dS /dV do not vanish but just move closer to zero, thus effectively condensing then all 
together in a single step at V — 0. The height of this step is not phase sensitive and is 
simply given by the expression (8) for the dc noise averaged over the period 2ir/Q. 

It is worth mentioning that our results for S u are quite general. Indeed, it is clear after 
what have been said that the singularities at V = nhQ/e are only due to the sharp edge of 
the Fermi distribution, and not related with any specific geometry assumed for the junction. 
Because of that the phenomenon should be displayed by any coherent conductor, provided 
that the main source of inelastic scattering is the ac potential. The reason is that an elastic 
scattering, if any, can smear the Fermi distribution of momenta but it will not affect the 
sharpnes of the step in the energies distribution, and our effect is sensitive only to the latter. 
The same remark applies to the oscillating dependence of the singularities on the amplitude 
of the ac signal. 




(8) 
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Let us briefly discuss a generalization of the system shown in Fig.l where the loop is not 
an ideal lead but a real metallic wire with disorder, i.e., instead of one scatterer there are 
now many of them uniformly distributed over the bulk of the wire. Most interesting is the 
case of a purely coherent conductor for which the energy relaxation time te and the phase 
breaking time are much longer than the flight time £/. (One can estimate £/ ~ Ti/E c , 
where E c is Thouless' energy hD/L 2 .) In such a system transport is described by channels 
of the scattering matrix with transmission coefficients T m assigned to each channel M. In 
the dc case the noise can be written || in terms of T m as So = — J2 m T m (l — T m )eV. In the 
presence of the alternating flux the extension of our formalism can be carried out easily and 
one obtains expressions similar to (6) and (7), with D 2 and D(l — D) replaced by S m T^ 
and Em^mfl — T m ) respectively. However, the limitations under which the result is valid, 
eV <C h/tf, ft <C 1/tf, are now slightly more stringent than for Exprs.(6),(7) because the 
flight time tf is longer. 

A more fundamental limitation to the general validity of our calculation is in the assump- 
tion that the flux threads only the phase coherent part of the conductor. It would certainly 
be of interest to better understand the opposite limit when the ac voltage increases smoothly 
over a distance much larger than the phase breaking length = ^Jt^/D. 

To summarize, we studied current and noise in a conductor driven by dc and ac and we ex- 
pressed them through time- dependent one particle scattering amplitudes. In the quasistatic 
limit of short time of flight through the conductor the current is given by the Landauer 
formula with time- dependent transmission coefficient, i.e., by a trivial generalization of the 
static case. The situation with the noise is quite different because of the two-particle in- 
terference. Spectral density of the noise depends on the scattering amplitudes in such a 
way that the phases do not drop out, and this leads to a non- stationary Aharonov-Bohm 
effect. Because of the way the Fermi statistics affects the two-particle interference the noise 
measured at T = is singular at u — ±eV/h + mQ, where m is any integer. To illustrate 
the phase sensitivity of the noise we consider a conducting metallic loop in which the ac 
signal is supplied by an oscillating magnetic flux. Because of the sensitivity to the phase 
of transmission amplitude the strengths of the singularities in the noise display oscillatory 
dependence on the amplitude of the ac flux given by squares of the Bessel functions. 
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FIGURES 

FIG. 1. Differential noise dSo/dV at T = given by Exprs.(5),(7) is plotted against V for three 
flux amplitudes: (1) <J> a = 5$ /4vr; (2) $ a = 7$ /2vr; (3) $ a = 23$ /47r. (Inset: Junction with 
leads bent in a loop through which alternating magnetic flux is applied.) 
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